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A LOCAL VERSION OF THE HUGHES MODEL FOR 
PEDESTRIAN FLOW 

JOSE A. CARRILLO*, STEPHAN MARTIN * 1 ", AND MARIE-THERESE WOLFRAM 1 


Abstract. Roger Hughes proposed a macroscopic model for pedestrian 
dynamics, in which individuals seek to minimize their travel time but 
try to avoid regions of high density. One of the basic assumptions is that 
the overall density of the crowd is known to every agent. In this paper 
we present a modification of the Hughes model to include local effects, 
namely limited vision, and a conviction towards decision making. The 
modified velocity field enables smooth turning and temporary waiting 
behavior. We discuss the modeling in the micro- and macroscopic set¬ 
ting as well as the efficient numerical simulation of either description. 
Finally we illustrate the model with various numerical experiments and 
evaluate the behavior with respect to the evacuation time and the overall 
performance. 


1. Introduction 

The mathematical modeling and simulation of pedestrian dynamics, such 
as large human crowds in public space or buildings, has become a topic of 
high practical relevance. The complex behavior of these large crowds poses 
significant challenges on the modeling, analytic and simulation level. These 
aspects initiated a lot of research in the mathematical community within 
the last years, which we briefly outline below. Mathematical modeling ap¬ 
proaches for pedestrian dynamics can be roughly grouped into the following 
categories: 

(1) Microscopic models such as the social force mo del [231 (22l 36j or cel¬ 
lular automata approaches [10]. 

(2) Fluid dynamic approachesfTT\ [3li. 3] and related macroscopic models, 
see for example the popular Hughes model [27., Ij El EJ 37] . 

(3) Kinetic models [5j [16] which uses ideas from gas kinetics to models 
interactions between individuals via so-called collisions. 
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(4) In optimal control [25] and mean field game approaches \29\ 118] pedes¬ 
trians act as rational individuals, which adjust their velocity optimal 
to a specific cost. 

(5) Multiscale models coupling between different scales to describe for 
example crowd leader dynamics mm- 

A detailed survey on crowd modelling can e.g. be found in Bellomo and 
DobgefBj. Several aspects are considered to be important in the mathemati¬ 
cal modeling to capture the complex behavior in a correct way. For example 
repulsive forces when getting too close to other individuals or obstacles play 
an important role in the dynamics. Another popular assumption is the fact 
that individuals act rationally and try to make the optimal decision based 
on their actual knowledge level. Partial knowledge of the overall pedes¬ 
trian density or the domain is another important factor which should be 
taken into account in the modelling. While these nonlocal effects can be 
implemented quite intuitively on the microscopic level, their translation for 
macroscopic models is not straightforward. Most macroscopic nonlocal mod¬ 
els are based on the continuity equation for the pedestrian density, where the 
nonlocal effects correspond to the deviation of the crowd from its preferred 
direction |llj 121113]. This deviation is determined by the average density 
felt by the pedestrians and modelled via a convolution operator acting on 
the velocity. The development of numerical schemes for conservation laws 
with nonlocal effects gained substantial interest in the last years. This was, 
among other factors, also initiated by the development of nonlocal models 
in traffic flow Els]. 

The original model of Hughes [27] describes fast exit and evacuation scenar¬ 
ios, where a group of people wants to leave a domain SlcK 2 with one or 
several exits/doors and/or obstacles as fast as possible. The driving force 
towards the exit is the gradient of a potential = 4>(x, t), x 6 fl, t > 0. This 
potential corresponds to the expected travel time to manoeuvrer through the 
present pedestrian density towards an exit. Hughes assumed that the global 
distribution of pedestrians is known to every individual, an assumption not 
generally satisfied in real world applications. 

In this paper we present a generalization of the classical Hughes model, 
which includes local vision via partial knowledge of the pedestrian density. 
We discuss the proper modeling setup, the implementation of suitable nu¬ 
merical schemes as well as their computational complexity. Furthermore we 
compare how the reduced perception of each pedestrian effects the overall 
"performance" of the crowd in evacuation scenarios. Inevitably, one expects 
the crowd to behave less efficient as less information is available. Quanti¬ 
fying how localised vision influences performance and decision making is a 
very interesting question in terms of collective behaviour. Surprisingly, it 
will turn out that evacuation times can even improve. The question we in¬ 
vestigate is therefore complementary to mean-field game approaches, where 
pedestrians anticipate future crowds states and hence are more capable than 
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the original Hughes’ model [29, T8] f9]. 


This paper is structured as follows: We start with a review on the modeling 
and analytic results of the classical Hughes model for pedestrian flow and 
its microscopic interpretation in section [2] In section [3] we present the local 
version of the Hughes model on the micro- and macroscopic level. Section [4] 
presents the numerical strategies for the microscopic and macroscopic model. 
We compare the behavior and performance of the models in Section [5] and 
conclude with a discussion of the proposed model in Section [6] 


2. Hughes’ model for pedestrian flow 

2.1. Original formulation and analytic results. Let us start by pre¬ 
senting the original modeling assumptions and the corresponding partial dif¬ 
ferential equation system of the Hughes model for pedestrian flow. Hughes 
considered an exit scenario, in which a crowd modelled by a macroscopic 
density p = p(x, t ) wants to leave a domain as fast as possible. The nonlin¬ 
ear PDE system for p and the potential cj) = t) on the domain Hcl j 

reads as: 


(la) 

(lb) 


^ -div(p/(p) 2 V</>) = 0 

IIWH = 


1 


f(p)' 


The first equation describes the evolution of p in time, driven by the gra¬ 
dient of (f> and weighted by a nonlinear mobility / = f(p). This mobility 
includes saturation effects, i.e. degenerate behaviour when approaching a 
given maximum density p max £ Possible choices are f{p) = p — p max 
or f(p) = (p — /0 max ) 2 amongst others. The former is inherited from the 
Lighthill-Whitham-Richards model for one-dimensional traffic flow [301 1 33] . 
The potential (j) corresponds to the weighted shortest distance to an exit in 
the following sense: Solving the eikonal equation (llbl) determines the op¬ 
timal path V0 minimising the expected travel time throughout the crowd 
towards an exit. This cost is measured as the inverse of f(p), hence the cost 
of walking through dense regions is high. Equation fllbl) is also a stationary 
Hamilton- Jacobi-Bellman equation, and the optimal path property of V0 can 
be rigorously derived m [ 24 ]. The fact that the potential cj) solely determines 
the direction of the flow can be easily seen as / 2 (p)V</> = f(p)\7c/)/\\ V0|| 
using 

Hughes model (jTJ) is supplemented with different boundary conditions for 
the walls and the exits. We assume that the boundary is divided into two 
parts: either impenetrable walls 9H wa u C or exits/doors <9fl ex it C 
with <9H wa ii n <9fi ex it = 0- Typical conditions for the density p in (JT]) are zero 
flux boundary conditions at <9fi wa u, which are either automatically satisfied 
as V</> • n = 0 or artificially enforced. The flux at 3H ex it has to be defined 
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according to the arriving density and our choices are discussed in Section [3] 
The boundary conditions of fllbl) are set as 4>(x,t) = 0 for all x G <9f2 ex it- 
There has been a lot of recent mathematical research on the classical Hughes 
model PHEBfflEO]. Up to the authors knowledge all analytic results are 
restricted to ID only, which is caused by the low regularity of the potential 
4>. This low regularity, i.e. (f> G C 0,1 , results in the formation of shocks 
and rarefaction waves in the conservation law. It is caused particularly by 
the existence of sonic points, which are hypersurfaces in space, where costs 
towards two or more exists coincide, and therefore V</> does not exist and 
the orientational field is discontinuous. In spatial dimension one the system 
can be reduced to the conservation law with a discontinuous flux function. 
In this case it is possible to solve the corresponding Riemann problem HZ3, 
which also serves as a basis for different numerical schemes GUIS]. 


2.2. Microscopic interpretation. Hughes motivated system © on the 
macroscopic level only. Recently Burger et al. |9j were able to give a micro¬ 
scopic interpretation of ©, which will serve as a basis for our local particle 
model. Microscopic models based on Hughes’ modelling assumptions are 
also used in the field of computer vision |36j. 

Let us consider N particles with position X 3 = X 3 (t) and velocity V 3 = 
V 3 (t), j = 1,... N. Then the empirical density p N = p N (t ) is given by 

1 N 

(2) p N (t) = -Y,Kx-X , (t)). 

3 =1 

Furthermore we introduce its smoothed approximation p^j = p^ (t), given 

by 


Pg(x,t) = (p N *g)(x,t) = - x 3 (t)), 

3 = 1 


where the function g = g(x ) corresponds to a sufficiently smooth positive 
kernel. The walking cost is given by the sum of a weighted kinetic energy 
and the exit time, defined as T ex n = sup{t > 0 | x G H}. Then the problem 
reads as: 


(3) 


C(X;p(t )) = min - 

V ry (x,v) 2 

subject to — = V(s) and £(0) = X(t). 
as 


^ H y(s) H 2 dc I V (XV) 
t W(£M),*)) 2 lex ' t[X ' Vh 


Hence the optimal trajectory is determined by ’freezing’ the empirical den¬ 
sity = p N (t ), in other words it corresponds to extrapolating the empirical 
density p N in time when looking for the optimal trajectory. 

Burger et al. were able to show that Hughes’ model can be formally derived 
from the optimality conditions of © and letting T —» 0 (corresponding to 
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the long-time behavior of the corresponding adjoint Hamilton-Jacobi equa¬ 
tion) . 

We will use this microscopic interpretation to propose a numerical approx¬ 
imation by a particle method in Section 4 of Hughes type models. In fact, 
(jlal) is seen as a continuity equation with velocity field v(x,t) = —/(p) 2 V^> 
driven by fllbl) . and thus particles in (j2|) are advected by the velocity field v, 
e.g. 

dXi 

— =v( X 3{t),t), j = 1, . . . N . 

3. A LOCALIZED SMOOTH HUGHES-TYPE MODEL FOR PEDESTRIAN FLOW 

The Hughes model (JT|) assumes that at any time t > 0 the global distribution 
of all other individuals p(x,t ) is known to every pedestrian. Therefore she 
chooses her optimal walking direction V</> in order to minimise its expected 
travel time/costs. Here, all walking costs are based on the current density, 
which means that pedestrians do not anticipate future dynamics of the crowd. 
Instead they are capable to react to changes in the global density ad-hoc as 
the path optimisation is repeated continuously in time. In a mean-field game 
type model, the capabilities of pedestrians would increase, as the planning 
decision of all agents can be correctly predicted into the future. We follow 
an opposite approach and reduce the capabilities of pedestrians, to obtain a 
more realistic model. 

The assumption of continuous and complete perception of global density 
information at current time is highly questionable in practical situations. 
Limited vision cones and restricted perception of global information comes 
through obstacles (walls, buildings), physical distance, visual orientation or 
the inability to see through a very dense crowd. Some effects of local vision 
on the behaviour of crowds are obvious: in an evacuation scenario with two 
exit corridors, which cannot be seen from each other, pedestrians caught in 
a jam in front of one exit will not be able to see whether the other exit is 
free or also jammed. 

These considerations motivated a new version of Hughes-type pedestrian 
dynamics based on localised perception of information, which we introduce 
in this section. The decision of each pedestrian is based on the perceived 
local density available in a limited domain, which can be e.g. interpreted as a 
vision cone. Furthermore a local interaction mechanism between individuals 
as well as a smoothening kernel on the velocity field (to prevent unrealistic 
high frequency oscillations in the direction of motion) are incorporated. We 
begin with the detailed introduction of the macroscopic model and discuss 
the microscopic analogue thereafter. 

3.1. Macroscopic equations. The starting point of our model is the as¬ 
sumption that pedestrians still perform the same path-optimisation selection 
as in the Hughes’ model, while the crowd state they act upon subjectively 
depends on their position and the amount of information they are able to 
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perceive. Let y € M 2 be an auxiliary variable and (j)(x, y ) : K 4 -A- R a 
parameterised potential, such that y i-a <f>{xQ,y ) denotes the cost potential 
calculated by pedestrians located at xq £ 0. To model space-dependent per¬ 
ception of information, suppose that for every x the domain 17 decomposes 
into a visible subdomain V x 3 x and a hidden or invisible part H x = 

We propose the following mechanism of restricted vision: If an area is visible 
its density is known and priced accordingly in the path optimisation. If how¬ 
ever an area is not visible, its density is thought to be a constant pu £ Mq", 
which we assume to be uniform among all pedestrians. Exemplarily, setting 
PH = 0 implies that pedestrians assume that not visible areas to them are 
empty, hence they will have a strong incentive to explore unseen parts of 
the domain. On the contrary, pedestrians will avoid invisible areas when 
PH ~ Pmaxi as they assume high costs. An eikonal equation in 17 is hence 
solved in the auxiliary variable y for every point x as 




which gives the potential <j> as function of two space variables. Note, that 
this notion of local perception differs from other recent work[8], where a local 
average of the density is used. Each pedestrian uses the cost potential at her 
own position for the decision process. Computing V y <j)(x,x) hence would, 
after normalisation, give a new orientational vector field to be used in the 
unchanged transport equation. We however argue that it makes sense to 
include a notion of conviction to the model, which has previously not been 
considered. In order to do so, (J4]) is solved for every single exit. This results 
in the computation of M potentials cj)k = (j>k(x,t), k = 1 which 

allows for cost comparison between exits, see Remark 3.2. 

In regions of high density, decisions on the walking direction towards any 
of k = 1,... ,M exits <91! E k cannot arbitrarily deviate between neighbours. 
If a pedestrians prefers to walk against the direction of a predominant local 
flow, collision or friction losses in the movement will occur. Especially on the 
macroscopic level, in which we take an aggregate perspective, an incentive 
to change the flow cannot arise from one point in space alone. As we do not 
model the granular level of individual collision or friction, we propose the 
following mechanism: 

(1) Each pedestrian carries an individual conviction strength u(x) mea¬ 
suring its preference of its chosen exit over all others. 

(2) There exists a local consensus process within the crowd, which results 
in the adjustment of the individual walking direction according to the 
predominant direction around them. 

Hence, pedestrians adjust their own direction in order to prevail the flow 
rather than obstructing it. This can be seen as either a cognitive decision 
rule or a forced physical restriction. For a compactly supported interaction 
kernel K : R 2 —> M, we define the final walking direction tp(x) at any point 
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x £ fi as 

(5) 


V = 


pU-kK. 
pkK ’ 


where the conviction ■u(x) is given as 


( 6 ) 


u = 


\7 y(f>kO P t 
| Vjy(/>fcopt || 


((j) k opt + 1 — (j) k opt), 


obtained by comparing the cost potentials 0*., k = 1,.. 
each of the exits: 


., M. associated to 


(7) k°P\x) = argmin fc ^(x, x), 

(8) ^°pt+i(x) = argmin^opt <j> k (x, x). 

Discontinuities in the velocity field due to the heterogeneity of decision 
making amongst pedestrians are hence partially compensated. To further 
smooth the model, we relax the strict restriction of ||V</?|| = 1 of Hughes’ 
model and replace the normalisation operator with a smooth approximation 
V : M 2 —> M 2 defined as: 


(9) V[x\ := < 


M 


||x|| > £, 

Sin ( 2 arctan(fcl) axctan(fc||x||)) ^ 0 < ||x|| < £, , 

0 x = 0, 


for some parameters k,£ > 0. We stress that this is not a technicality, as we 
here allow pedestrians to stop when being undecided. This is highly desirable 
from the modelling point of view, though on the other hand the modulus of 
the flux now is not a function of density alone, as one can see below. 

Next we discuss the boundary conditions for the eikonal equations. Since we 
treat each exit separately, we set = 0 i R the computation of <j) k . No 

boundary conditions are imposed on the rest of the boundary 3H wa n. 
Near-wall and near-obstacle effects have a strong influence on the dynam¬ 
ics on constrained macroscopic evolutions. We propose that pedestrians 
take into account walls and obstacles in their computation of optimal paths. 
Hence it is natural to include these effects as an additional fixed cost W (x) 
on the right-hand side of the eikonal equation ((4]) . We introduce a smooth 
layer profile Xw( x ) £ [0,1], which identifies areas close to walls but smoothly 
vanishes elsewhere and around exits to allow outflow. A typical choice of Xw 
is illustrated in Figure [T] For the sake of simplicity, we set 


( 10 ) 


W(x) = 

/ (Pmax Cj 


hence areas close to walls are penalised similar to high density areas. 

Finally all terms are coupled to the continuity equation with velocity held 
v(x,t) = —f(p)V[V(p\ as in the original model. At exits, we prescribe a 
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0 0 



Figure 1. Illustration of layer profile Xw : wall and obstacle 
repulsion are embedded to the model with a fixed cost W 
defined in terms of the layer profile Xw( x ) € [0,1], see (1101) . 
which indicates proximity to walls less than a width w, but 
vanishes away from obstacles and near exits to allow a proper 
outflow of pedestrians. Xw equals one at walls and is the linear 
function given by vertex values within triangles Tf, T 2 near 
exits. 


maximum outflow, given by v(£,t) = — f(p)n for all £ G cftl e xit- Taking all 
these considerations into account the full macroscopic model reads as: 


(lla) 

(lib) 

(11c) 

(lid) 

(lie) 

(Ilf) 

(Hg) 

(Hh) 


d t p(x, t) + V x - ( ~f(p(x , t))V[V<p(x , t)\p(x, t)) = 0 

(pu * IC)(x, t) 


<p(x,t) = 


(p*K,)(x,t) 


u(x,t ) = T 7 ^Z~T ~~~~ X ' ^ (0fc° pt+1 (x, x , t) - (j) kop t(x,x,t)) 

fc opt (x, t) = argmin fc cj) k (x, x, t ) 
fc° pt+ i(a;, t) = argmin fc ^ fe o P t 0fc(x, x, t ) 

1 7ife) y e Hx 

s-t.: = 0,fc = l,...,M,Vi 

s.t V[p{^, t)\=n, V£ G dQ exit , s.t.: p(x, 0) = p 0 (®)- 


l|Vj / 0 fc (a;,y,t)|| = 


We conclude the section with remarks on specific modeling assumptions. 


Rem.ark (vision cones). We have set aside formal statements regarding as¬ 
sumptions on the visible set V x , but clearly we think of at least regular, 
connected and closed sets. A necessary condition is 

(12) x G V° = int 14, 

which implies that every pedestrian perceives some information from all di¬ 
rections. This restriction rules out e.g. angular vision cones (see Degond 
et.al.[;15j) where pedestrians do not see what is happening behind them. 
In our model, f|T2l) is necessary to exclude unrealistic situations where the 
chosen walking direction points outside the visible area. The inclusion of 
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angular-dependent vision cones is certainly possible, but would imply a 
velocity-dependency and lead towards a second-order macroscopic model. 

Remark (conviction term). The introduction of the conviction term u{x) re¬ 
quires the computation of exit costs <f>k via the eikonal equation for individual 
exits, which appears to be a significant complication of the model. However, 
it is worth noting that the mechanism is almost identical to the original 
model. In equation fllbl) the costs of walking towards any of the K exits are 
compared, but only the minimal costs are used. Here, we simply store more 
information. This connection is also illustrated by looking at the numeri¬ 
cal schemes for solving the eikonal equation (see also Section ?): If a Fast 
Sweeping Method is used in e.g. a corridor with two exists, this essentially 
corresponds to solving for each exit separately if the minimization step is left 
out. If a Fast Marching Method is used, the conviction is directly related to 
the sequence in which vertices are promoted, with the least convinced vertex 
being assigned a cost the latest. 

Remark (waiting behavior). The relaxation ||'P[cc]|| < 1 implies that the 
modulus of the flux can be less than f{p)p when pedestrians are undecided. 
This makes a rigorous analysis of the model equations a difficult task, which 
is not tackled in this work. The benefit of our formulation is that the problem 
of discontinuous velocity fields at sonic points has disappeared. Pedestrians 
at those hyper surfaces will not move unless the sonic points move. 


3.2. The one-dimensional case. Consider a one-dimensional corridor Q = 
[0,1] with two exits and the uniform radial vision cone V x = [x — L/2,x + 
L/2] flO of length L. Exit costs towards the left and right exit are computed 
at y G 14 as 


4>l(x, y, t) = 


1 


-dz + 


J Kph) 

z<y,zdH x z<y,z£V 3 


/ f(p(z, 


t)) 


(f>R{x,y,t)= J 


l 


-dz + 


/ 


1 


Kph) 

z>y,zeH x z>y,z£V r a 

u(x, t ) = 4> l (x, x, t) - 4>r(x, X, t) 


f(p{z,t )) 


dz, 


dz, 


as illustrated in Fig. |2] The cost potential (j) is two-dimensional and d y (j)(x , y) 
gives the preferred walking direction that a pedestrian located at x seeing 
14 assigns to y G [0,1]. The walking directions chosen prior to the consen¬ 
sus process are hence given as d y (j)(x,x ) along the diagonal of [0, l] 2 . For 
every fixed x G [0,1], there is a unique sonic point z(x), where 4>l(x, z(x)) = 
4>r(x, z{x)) and d y (j)(x, z{x)) does not exist. As illustrated in Fig. [3] the indi¬ 
vidually preferred walking directions can switch multiple times between both 
exits, depending on the current density and the vision cones. At switching 
points, the preferred directions can point outwards (separation) or inwards 
(collision) and only the weighted interaction process (IllbMllcI) generates a 
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Figure 2. Illustration of path optimisation mechanism in 
ID: A pedestrian located at xq computes and compares the 
cost potential cj)i, of left vs. right exit in a corridor [0,1]. 

Next to its own negligible density, the present crowd consists 
of three blocks. Outside the vision cone V x , the evacuation 
costs grow linearly at constant rate, as the local density is 
unknown. Within V x , the slope of the cost potential increases 
with the pedestrian density. Preference is then given towards 
the exit with lower estimated cost. The conviction towards 
this decision is given as the cost benefit |n(xo)| = \4>l — 4>r\- 

smooth velocity profile. In the Hughes’ model, all vision cones are identical 
and there is a single separation point. 

3.3. Microscopic interpretation. We conclude this section by briefly com¬ 
menting on the modelling of local vision at the microscopic level. The mi¬ 
croscopic modification is straightforward and uses the same ideas as at the 
macroscopic level. It corresponds to updating the position X = X ( t ) accord¬ 
ing to a potential which depends on local information only. Its calculation is 
based on the same equations as in the macroscopic model f lllfl) but using the 
smoothed empirical density p^ instead of p. The position update is based on 
equations flllbll - flllcl) . Hence individuals choose the path towards the exit 
with the lowest cost, but weigh their decision according to the predominant 
direction chosen around them. For further details on the implementation we 
refer to Algorithm [2] presented in Section |4] 

3.4. Analysis of the domains of dependence. In this subsection, we will 
discuss some mathematical properties of the solutions of the eikonal equa¬ 
tions flllfl) . From the construction of the model, the potential cj)(x,y,t) has 
to be computed for every x £ 0 on the entire domain H, which counterbal¬ 
ances the idea of locality and increases the computational cost considerably. 
We show here, that the computation of the potential can actually be reduced 
to a subset of H, called the effective domain of dependence, for every x. Only 
this subset, which contains V x , is considered in the individual local planning 
problem and corresponds to the reduction of the computational cost. 





















A LOCAL VERSION OF THE HUGHES MODEL FOR PEDESTRIAN FLOW 


11 


1 



1 


X 


0 


0 


0 


y 


1 


0 


y 


i 


(a) Localized model 


(b) Hughes’ model 


FIGURE 3. Illustration of turning decisions of a ID popu¬ 
lation for a given p(x,t ): For every point x, we show the 
individual sonic point z(x), where the costs z(x)) = 

(f>n(x,z(x)) coincide. The preferred walking direction for a 
pedestrian at x is found along the diagonal (x,x). If the 
sonic point is to the left, the pedestrians aim to walk towards 
the right and vice-versa. No direction is preferred and the 
conviction is zero if the curve of sonic points intersects the 
diagonal, (a) Local vision cones: The preferred direction al¬ 
ters and creates multiple points of separation and collision. 

The resulting velocity field is obtained by the smoothening 
interaction process (IllbD - flllcI) . (b) Hughes’ model: All vi¬ 
sion cones coincide, hence there is one identical sonic point 
common to all pedestrians. 

The following proofs rely crucially on the optimal path property of the 
characteristics associated to the eikonal equation flllfl) . We reca 11 1231 S] that 
by Fermat’s principle the characteristic paths associated to (j)(x,y,t), given 
by the solution of: 


(13) 7*,t0) C H : 7 ( 0 ) = z, 7 (s) = - V0(x, 7(5), t) for all s > 0 , 


are the optimal paths for the cost defined as 



y e 14 
y e H x 


Moreover, the potential is the value function for that cost. Hence it is de¬ 
creasing along these paths and satisfies the optimality condition 



0 < a < b 


J a 


being zero at its corresponding exit <9f2 e xit- Furthermore, the curves 7 ^ t are 
the optimal paths to achieve the exit, i.e., they verify the following global 
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optimality condition 

rT z rTz 

(15) +<j>(x,z,t)= c('y% t (s),t)ds< c(j(s),t) ds, 

Jo Jo 

for all 7 curves joining z to any point in the exit dfl ex a, where T z is the 
optimal time to achieve the exit for the point z € fl and T z is the time to 
achieve the exit for the path 7. 

Lemma 1. Consider any fixed V x C Ll and that f(p) >0, 0 < p < /0 max - Let 
<f>H be the global solution of the eikonal equation ||V</>ij|| = ^-//(ph), <t>H{x) = 
0 on Sriexit- Define the minimum of fin in V x as 

rriH := min 4>h(z) 

and the corresponding superlevel set of fin as 

Mh := {x € : 4>h{x) > m^}. 

Then the problem of computing the local potential </>(x, y, t ) =: 4>{y) out of 
(lllfl) on 17 reduces to the following problem on M^\ 

'\\Vyky)\\ = J^+ W (y) m 14 

< \\Vy(i>(y)\\ = in M h \V x 

M:y) = m H on dM H \dn wa , u (b.c. ) 

Proof. If an exit is visible then rri h = 0, Mh = Ll and the assertion is trivial. 
If no exit is visible then by construction V x C Mh and 4>h = wh > 0 
on OMh■ As the walking costs are always positive, c(y,t ) > 0, we get 
4>(x, y, t) > mn for all y E int Mh- On the other hand, any point z E Q\Mh 
satisfies cf>(x,z,t) < mn and hence 7 x t( s ) does not intersect V x , otherwise 
the cost should be larger at a middle point than initially contradicting the 
optimality of the path 7 xt (s) in ( I14[) . Hence OMh is the maximal level 
set consisting of points whose optimal paths do not cross V x , and therefore, 
4>{x, z, t ) can be computed from flllfl) with constant right-hand side outside 

M H - □ 

Definition 1. Consider a fixed visibility area V x . For a z E H, denote the 
default optimal path 7^ as the parameterised curve associated to a gradient 
walk along fin starting in z. that is 

T h (s) C H : 7(0) = z,j(s) = -V0h( 7(«)) Vs > 0. 

Next, dehne the characteristics’ shadow V& as the set of all points, whose 
default optimal paths crosses the visibility area, hence 

V# := {z E : 7^ C int V x 0}. 

Note, that V V C Mh since any default optimal path outside of Mh cannot 
intersect with V x as proven in the previous lemma. 
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Lemma 2. Consider any fixed V x C. Ll and assume that f(p) > 0 is increasing 
in 0 < p < /Jmax, with pu = 0, then the problem of computing the local 
potential <j)(y) out of (|4|) further reduces to the following problem on V& 

f \\Vyky)\\ = 7 ^ + w (y) in 14 
< II VyHy)\\ = j^ mV#\V x . 

^ 0 = 4>h on dV# 


Proof. For any point z whose default optimal path 7^ that does not intersect 
with V, the claim is that (j)(z) = 4 >h( z ) due the monotonicity of the cost 
function, i.e., 

————- + W(y) > —— . 

f(p(y,t )) /( 0 ) 

To prove this, let us denote by Tf 1 the optimal time to get to the exit for 
the default optimal path 7 fj. 

We first take 7* t (s) as a candidate path in the global optimality condition 
m for the eikonal equation with right hand side ch = jjaj- Being 7 xt (s) a 
path joining z to a point in the exit and 7 h(s) the optimal one, we conclude 

rT z 


<I>h(z) = T^c h < T z c h < I c(7j ,(s),f) ds = f>(z). 


f 

Jo 


Now, we take 7 h( z ) as a candidate path in the global optimality condition 
(1151) for the eikonal equation with right hand side c(y,t). It is an admissible 
path as it connects z to a point at the exit and the cost along its path 
coincides with ch = c(7|j-(s),f) for all s £ [0 ,T%] since the path does not 
cross V. Then, we get 


<t>(z) < 


f 


c(l z H (s),t) ds = T^ch = <I>h(z) 


leading to the stated result. 


□ 


We illustrate Lemmata [Tj and |2] in Figure [4] It can be seen, that the 
reduction of the computational domain from Mh to V& can be significant, 
as the size of Mh depends on the closeness of V to the nearest exit, not 
on the size of V. For the exemplary geometry of Figure [H the boundary of 
V# coincides with the sonic points of 0//, but this is not true in general. 
Furthermore, it is easy to see why the computational domain cannot be 
reduced further. Suppose that p(-,t ) is spatially homogeneous, then —V0 in 
V#\V points to the left exit as in the eikonal case. On the other side, one 
can choose a situation with a large density at the left boundary of V that 
leads to right-pointing —V0 in V#\V. 


4. Computational methods 

In this section we present a microscopic and a macroscopic numerical 
solver to simulate the classic and the local version of the Hughes model. The 
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Figure 4. Illustration of the domains V,Mu and V # for 
computation of the visibility area potential (j) for the case 
of a corridor with two opposing exits: The problem on 17 
generally reduces to a HJ-equation on Mu ■ as by construction 
$ coincides with (f>u outside of Mu (Lemma (Jl]), — V<^># solid 
arrows). If pu = 0, any default optimal path of (j)u that 
does not intersect V remains optimal, as indicated by dotted 
arrows, and the problem reduces to V# (Lemma 0. 


proposed methods have been implemented on regular and triangular meshes 
in 2D to allow for flexible discretizations of polygonal domains with one or 
several obstacles. 

For the macroscopic system flTTl) we use the following explicit iterative 
algorithm: 


Algorithm 1 Macroscopic version of the localised model 

Initialisation: 

• A discretisation Cl = (V, £, T) of Cl consisting of vertices, edges and 
cells. 

• An initial density po given on T, such that f T p(x, 0 )dx = 
Po(T) VT6T. 

• A list of exits, a list of boundary edges per exit and |V| subsets of V 
containing the vision cones defined in terms of vertices. 

(1) Compute the cost potential (pk for all exits out of the current density 
p by solving (lllfl) along the vertices for every v £ V. 

(2) Determine the cell values of (f>k and V0& by an averaging / finite 
difference approximation of the values at neighbouring vertices, e.g. 

= \{vedT}\ YlvedT&kiu), and obtain u(T ) here from. 

(3) Compute a numerical convolution of u with 1C, which gives (p on the 
cells. 

(4) Update the density with a cell-based Finite Volume Method using 
the velocity held — f{p)'P[0\ and a suitably chosen time step. 


The discretisation is either a regular grid or an unstructured regular tri¬ 
angular mesh to allow more complex geometries. For solving the eikonal 
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equations, one can chose between Fast Sweeping Methods|381 i32] and Fast 
Marching Methods[28j [34]. The former is based on a Gauss-Seidel iteration, 
which updates the solution by passing through the computational domain 
in alternate pre-defined sweeping directions. A rectangular grid provides a 
natural ordering of all grid points. This ordering does not exist on an un¬ 
structured grid and is replaced by a general ordering strategy by introducing 
reference points, which is done once. Then the solution at each node is con¬ 
secutively updated by running through the ordered lists. Marching methods 
update vertices in a monotone increasing order, where in every iteration a 
list of candidate values is available by finite difference approximation from 
previously approved values. The smallest value all of candidate values is 
then promoted and assigned to its vertex. 

As a Finite Volume Method we use the first-order monotone FORCE 
scheme[55. [19] . Some postprocessing between the steps of Algorithm [T] is re¬ 
quired: outward-pointing components of V</>fc are removed along the bound¬ 
ary, suitable values of V0& are ensured at corners of D, and the max outflow 
condition flllhl) is enforced at cells neighbouring exit edges. 

The analogous algorithm used for the numerical simulation of the micro¬ 
scopic model is: 


Algorithm 2 Microscopic version of the localised model 

Let us consider a system of N particles, which are initially located at 
positions X 3 (0) = Xq . In every time step t l = iAt we update the particle 

position as follows: 

(1) Determine the empirical density at time t l : 

1 N 

Pg (X, = g ( X ~ X ^ (**))> 

i-i 

where g denotes a Gaussian. 

(2) Solve the eikonal equation to determine the weighted distance to each 
exit (j) k = </> fc (x,f*), k = 1 ,.. . M: 

+ W(y) if y € V x 
otherwise. 

(16b) <j) k (x,t l ) = 0 

(3) Update the position of each particle X J via: 

(17) X3(?) = • V^,f))), 

where (p(x,t l ) is determined by (jllbl) . 


(16a) 


||V^(x,y,f)|| = 


/ 


All _ ) f{pq{y V)) 


l 


Kph) 
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5. Results 


In this section we illustrate the dynamics of the localised model for crowd 
dynamics with examples in one and two dimensions. In all simulations we 
consider an evacuation scenario of a corridor, where a given initial distribu¬ 
tion of people tries to leave the rectangular domain through either one of 
the two exits as fast as possible. We compare the evacuation time, i.e. the 
time at which all individuals have left the domain, with respect to different 
parameters, e.g. vision cones. In the case of a global vision cone we obtain 
Hughes’ type dynamics. As a flux law, we chose the LWR function 

(18) f(p) = pO-~p), 

setting /9 max = 1 throughout this section. 


5.1. ID corridor - macroscopic model. In our first example the domain 
H corresponds to the unit interval H = [0,1] with two exits located at either 
end, i.e. at x = 0 and x = 1. We consider an evacuation scenario in which 
two groups, one of them being densly packed, want to leave through either 
one of the exits: 


(19) 


Po(x) 


' 0.85 0 > a; > 0.3 
< 0 0.3 < x < 0.6 , 

0.25 0.6 > x > 1 


and we set the width of the vision cone to L = 0.75. The resulting dynam¬ 
ics are illustrated at 4 time steps in Fig. [5] Within the left block, some 
pedestrians decide to walk towards the right exit, as they are aware of the 
high density on their left and account for a higher walking cost compared to 
the relatively empty right hand side. After the separation the right-moving 
part evolves as a rarefaction wave, as known from the LWR model. As the 
distance between the wave and the left-moving shock grows, the effects of 
the local vision cone become apparent. At some point pedestrians moving 
to the right do not see the high density at the left exit anymore and start to 
turn around. Therefore the rarefaction wave splits again - one part contin¬ 
ues while the other one turns around and moves back to the left exit. The 
turn-around occurs in several stages: 

(1) A new sonic point arises, where pedestrians are undecided between 
both exits. The walking direction is unchanged as the local consensus 
process (11 lbh prevents an immediate switching. 

(2) When a critical mass of density and conviction opting for walking 
to the left, the velocity after consensus switches continuously and 
passing through zero. This creates a temporary collision point, as 
there a still pedestrians to the left of the sonic point which walk 
towards the right. 

(3) The density at the collision point increases, which causes pedestrians 
to the left of the collision point to turn around too, as a higher density 
is in their way, as it can bee seen in Fig. |5](c). 
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0 
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(a) t 0 = 0 



(c) t 2 = 0.71 



(b) ti = 0.31 


/I 


l 

L 

l 

j / 

1_A 


1 1 


0 0.25 0.50 0.75 1.00 


(d) t 3 = 1.29 


FIGURE 5. Exemplary evolution of the ID model showcasing 
a turnaround behavior due to localised perception of infor¬ 
mation [density p solid (—), speed v = —f(p)V[Vip] dashed 

(-) and directional conviction — cf>L dotted (•■)]. (a) 

Piece-wise constant initial density. Part of the left crowd 
initially decides to move right in order to avoid the high- 
density jam. (b) The separated block moves to the right in 
a rarefaction-wave manner, (c) The wave is again separated 
as the high-density jam gets out of sight for centrally located 
pedestrians, who hence prefer the left exit and turn, (d) The 
turnaround is complete and remaining pedestrians will exit 
on the left. 


(4) Finally, all pedestrians to the left of the initial sonic point have turned 
and walk towards the left ( Fig. [5][d)). 

This new behavioural pattern is entirely consistent with the idea of con¬ 
stant re-evaluation of the optimal path based on restricted information and 
cannot be observed in the original Hughes’ model. We note that without 
the smoothening properties of the model around points of equal costs one 
obtains strong oscillations in the turning behavior, which causes severe nu¬ 
merical problems. The exact parameters of the simulation can be found in 
Appendix IA.1I 
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(g) Exit percentage 


FIGURE 6. Dynamics of the microscopic model in a two- 
dimensional corridor with two exits at the right and left. We 
observe a similar behavior as in the ID simulation in 15.11 - in¬ 
dividuals (visualized by red triangles) initially decide to move 
to the more distant exit, but after the congestion at the re¬ 
solves in time, they turn around and take the closer exit, (a) 
- (f): particle solution for different times and L = 0.25, (g): 
exit percentage over time for different values of L. 


5.2. 2D corridor - microscopic model. We illustrate the dynamics of the 
microscopic model in a two-dimensional symmetric corridor Q = [0,1] X [0, 
with exits at the left and right side, i.e. x = 0 and x = 1. The ID case 
of section 15.11 can be interpreted as a projection of this two-dimensional 
geometry. We consider the same initial distribution of individuals, i.e. the 
positions of all 500 particles are distributed according to the initial pedestrian 
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density fH9l) . For L = 0.25, Figure |6][a)-(f) nicely illustrates a similar turn¬ 
around behavior as in the ID macroscopic simulations. At the beginning the 
group close to the left exit splits, one part exits through the left exit the 
other one moves towards the more distant right exit. As the density close to 
the left exit decreases in time, the group moving towards the more distant 
exit splits again, i.e. parts of the group turn around and move back again. 
We marked all individuals, which initially moved towards the right but then 
turn around, with red triangles. Furthermore, Figure E^g) shows the change 
of the evacuation performance for different sizes of the local vision cones 
L. Here we plot the percentage of the total initial mass outside the domain 
versus time. Decreasing L and hence the perceived information, we observe 
that the overall evacuation performance first is merely diminished, and only 
begins to drop significantly after a certain threshold. The evacuation time 
will approach the uninformed eikonal case L = 0, which is not shown. All 
parameters can be found in Appendix I A. 2 1 


5.3. 2D non-symmetric corridor - macroscopic model. Now we turn 
to the macroscopic model in two dimensions. Again we consider the corridor 
fl = [0,1] X [0, ^], the exits however form only a part of the left and right 
edges, hence we obtain a fully two-dimensional dynamics where boundary 
conditions matter. The left exit is located between (0,0) and (0, 0.1) and the 
right exit is the segment connecting (1, and (1,0.4). The initial density 
Figure [TJa) is given as a low density group of pedestrians on the left and a 
high density group on the right 


( 20 ) 


Po{x,y) 


'0.1 0.05 < x < 0.3, 0 < y < 0.25, 

< 0.95 0.6 < x < 0.95, 

0 otherwise. 


We first study the case of global vision L = oo V x = H in Figure [Tj In 
(b), the low density group turns towards left and is quickly vacated. The high 
density group on the other side splits along a curve of sonic points. Pedes¬ 
trians turning to the right cause a jam in front of the right exit, whereas 
left-turning pedestrians occupy the corridor in a rarefaction-type manner 
inherited by the physical flux law (1181) . Upon arrival at the left exit, pedes¬ 
trians pile up and form a new jam (c). Hence, a fraction of the density turns 
around again and heads for the right exit (d, around (0.5, .25)), having to 
cross most of the corridor again (e). However, most of the pedestrians are 
committed to the left exit and do not turn, because the severeness of the 
left jam does not compensate their expected travel time, and the left exit is 
vacated later than the right exit (f). 

Compared to the classical Hughes model, the relaxation term Q and 
the conviction-based interaction (lllbll - flllcl) allow for a smooth turning be¬ 
haviour. The wall-repulsion (HOI) causes a density gap, which is different 
to zero-flux conditions as in [26], but prevents any spurious effects of the 
boundary flux. 
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(e) Time t = 2.1 (f) Time t = 2.75 


Figure 7. Two-dimensional macroscopic dynamics : We 
simulate model fllll) with global visual perception V x = Q. 
Two groups of pedestrians are initially placed in a corridor 
with a lower left and an upper right exit. The high density 
group separates according to the path optimization mecha¬ 
nism, as illustrated in several time snapshots of the density 
in (a) - (f). The right exit is vacated before the left exit and 
the final evacuation time (not shown here) is ss 3.4. 


It is clear that even the planning algorithm incorporated into the classic 
Hughes’ model does not lead to an optimal evacuation. One reason for sub- 
optimality in Figure [7] is that pedestrians have to keep in motion constantly 
but cannot predict the occurrence of future jams. Hence, pedestrians are 
likely to walk towards an exit that will be blocked in the future, as seen in 
the example. 

In Figure |8] we study the same initial configurations with localised percep¬ 
tion and a radial vision cone of diameter L = 0.75. The initial separation 
phase (a) is similar to Figure [7] As pedestrians move from the right to the 
left, the right jam gets out of sight and its influence diminishes. At the same 
time, the density on the left becomes visible. At a certain point a balance is 
achieved and pedestrians locally accumulate around an area of equal walking 
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Figure 8. Two-dimensional macroscopic dynamics : We 
simulate model flTTl) with local visual perception L = 0.75 
and the same initial configuration as Figure [7] In the time 
snapshots of the density, we observe again an initial separa¬ 
tion of the high density group (a). A waiting phenomena of 
up to two local groups is clearly observed (b-d). Crucially, 
waiting pedestrians are able to choose their exit later (d-e), 
which leads to a rather simultaneously clearing of both jams 
(f), as opposed to Figure [TJ The final evacuation time (not 
shown here) is ~ 3.025. 


costs, where in this case they are able to stop (b-c). Hence, we observe a 
waiting behavior which cannot be observed in classical Hughes’ type models. 
Looking from (c) to (d), a high density jam forms at the left exit, which 
causes part of the left-walking pedestrians to turn right after enough con¬ 
viction is gathered. Together with some outflow of the first waiting group, 
a second waiting group is formed (d). Pedestrians in a waiting group choose 
to move if one direction becomes favorable. As both jams at the exits reduce 
at the same rate, the left waiting group walks to the left and vice versa (e). 
Finally, the waiting groups dissolve and the exits get vacated at a rather 
similar time (f), and the evacuation time improves compared to Figure [7] 
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FIGURE 9. Evacuation time of the macroscopic model for 
varying vision cone diameters L: The performance can im¬ 
prove with limited vision. L = 0 corresponds to the eikonal 
case whereas L > 2.5 implies unlimited vision in the given 
corridor. 

The fact that evacuation performance can improve under limited percep¬ 
tion of information is surprising at first glance. Our simulations give an good 
explanation for the phenomena: As pedestrians show a waiting behaviour, 
they are less likely to be trapped in the jam arising at the left exit. In 
fact, the waiting is made possible by the combined effect of multiple sonic 
points due to local vision and the smoothed turning mechanism. Naturally, 
this is not generally the case and cannot be a-priori answered. For small 
vision lengths L, the dynamics will converge to the velocity field given by 
the eikonal equation, which is our initial configuration will exit almost all 
pedestrians using the right exit and perform poorly. In Figure [9l we study 
the evacuation time of 99% of the initial mass as a function of the diameter 
L. We unexpectedly find in this case wo optimal values of L for which the 
evacuation time is minimal. The classical Hughes’ evacuation time (L large) 
is always less or equal than the eikonal case (L = 0), however there is no 
way to generally argue that there will always be a minimum in between. 

6. Conclusion 

In this work we introduced a localized smooth variant of Hughes’s model 
for pedestrian crowd dynamics. We regularised the original model, com¬ 
posed by an eikonal equation and a continuity equation. First by a local 
interaction term, which intermediates individual pointwise path optimisa¬ 
tion towards conviction-weighted walking directions. Secondly, we allowed 
pedestrians to stop, if they are undecided, using a smooth approximation of 
the normalization condition. Most importantly, we restrict the information 
on the global density each pedestrian can use for her planning algorithm to 





A LOCAL VERSION OF THE HUGHES MODEL FOR PEDESTRIAN FLOW 


23 


a local surrounding area. This is a very realistic assumption for large crowds 
that has not been considered in the literature so far. We presented both 
a microscopic and a macroscopic version, and illustrated the model compo¬ 
nents in the one-dimensional case. In terms of analytical results, a rigorous 
theory for these kind of equations in multiple dimensions is currently out of 
reach to the best of our knowledge. However, we were able to identify some 
qualitative properties of the dependence of the optimal path on the vision 
cone that allow for a reduction of complexity. 

The numerical approximation of the model on both levels has been dis¬ 
cussed and utilizes several techniques including sweeping and marching meth¬ 
ods, particle approximations and finite volume schemes. Though the numer¬ 
ical costs of computing a solution have increased due to the inhomogeneity 
of vision cones, we observe new effects and phenomena in the model based 
on our simulations. First, local groups of pedestrians are able to change re¬ 
peatedly their walking direction towards an exit. This ’multiple turn-around 
behaviour’ can explained by the multiple sonic points of the estimated walk¬ 
ing costs, which by construction cannot occur in the classical case. We stress 
that the smoothening and conviction terms are crucial to allow a swift turn¬ 
ing behavior, which is not trivial to model in first order equations. Second, 
the model replicates a waiting behaviour in case of undecided pedestrians, 
i.e. in areas where locally estimated walking costs towards different exits 
are equal. Surprisingly, we found that this waiting phenomena induced by 
localized information can improve the overall evacuation performance of the 
crowd. In our numerical example we observed two local minima when varying 
the vision cone diameter. 

To conclude, we have demonstrated that local vision effects can be im¬ 
plemented into first order models for crowd dynamics. This leads to new 
unforeseen phenomena and complex behavior, whose partial understanding 
via qualitative properties is important for the applicability of such equations 
to social-economic problems. On the other hand, this work illustrates the 
limitations to first order models such as Hughes’, where planning decisions 
are instantaneously updated and no social or cognitive memory is taken into 
account. From our point of view Hughes’ type equations constitute an im¬ 
portant building block for crowd models and a mathematically important 
object of study, but it cannot be expected to be fully realistic. 
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Appendix A. Simulation parameters 

A.l. Id macro parameters in Section 15.11 The macroscopic simulation 
in ID was implemented in MATLAB. The domain [0,1] was uniformly dis¬ 
cretized with Ax = 10' 4 . The time step was set to At = 5 • 10~ 5 . The vision 
cone was defined as V x = [x — L/2,x + L/2] D [0,1] with L = 0.75. The 
radial interaction kernel /C was chosen as the indicator function on the in¬ 
terval [0,0.05]. The smoothed projection operator was chosen as in (j9|) with 
t = 0.05 and k = 25. The wall repulsion W (x) was neglected in ID. Absorb¬ 
ing boundary conditions were applied at both exits. The cost function was 
numerically bounded at c(p) < 10 4 . 

A.2. 2d micro parameters in Section 15.21 The microscopic simulations 
are implemented using the software package Netgen/NgSolve. The domain 
was discretized in 1438 triangles, the time steps were set to At = 10“ 2 , 
the final time to T = 1.5. At time t = 0 we distributed the 500 particles 
according to the initial datum po used in subsection 15.11 The empirical 
density p^ was calculated using Gaussians with variance a = 0.05 for the 
smooth approximation g. The width of the local vision cone was set to 
L = 0.25 in Figure [6j a)-(e) and to L = 0.75, 0.5 and L = 0.25 in Figure [6ff ). 

A. 3. 2d macro parameters in Section 15.31 The macroscopic simula¬ 
tion in 2D was implemented in MATLAB. The domain = [0,1] x [0, /] was 
uniformly discretized with Ax = Ay = 10' 3 . The time step was set to 
At = 5 • 10' 3 . A Fast Sweeping Method was used to solve the eikonal equa¬ 
tions. The vision cone was defined as V x = {y : \\y — x ||2 < 77 } n Cl with 
varying diameter L. The radial interaction kernel 1C was chosen as 



with b = 0.05. The smoothed projection operator was chosen as in (0 with 
t = 0.05 and k = 25. The wall repulsion W(x) was defined with the width of 
the boundary layer function x.w set to w = 0.025. The wall costs W{x) were 
numerically bounded at W (x) < c(0.975). The cost function was numerically 
bounded at c(p) < 10 3 . The numerical accuracy for vanishing density was 
set to 10''. 
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